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In this paper, we explore the implications of a two-point discretization of an extra-dimension in 
a five-dimensional quantum setup. We adopt a pragmatic attitude by considering the dynamics of 
spin-half particles through the simplest possible extension of the existing Dirac and Pauli equations. 
It is shown that the benefit of this approach is to predict new physical phenomena while maintaining 
the number of constitutive hypothesis at minimum. As the most striking feature of the model, 
■ we demonstrate the possibility of fermionic matter oscillations between the two four-dimensional 

sections and hyper-fast displacements in case of asymmetric warping (without conflicting special 
relativity). This result, similar to previous reported ones in brane- world theories, is completely 
£SJ ■ original as it is derived by using quantum mechanics only without recourse to general relativity and 

bulk geodesies calculation. The model allows causal contact between normally disconnected regions. 
If it proves to be physically founded, its practical aspects could have deep implications for the search 
of extra-dimensions. 
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I. INTRODUCTION 



The idea that our observable universe could be a part of a more extended N-dimensional spacetime (N>4) has a long 
tradition. It traces back to the seminal work of Kaluza in 1921 who extended general relativity in five dimensions in 
order to treat electromagnetism and gravitation on an equal footing [1] . Unfortunately the model and its subsequent 
extensions reveal unsuccessful in their aim of describing the physical reality such that the multidimensional approach 
was abandoned for a while. 

Since that time, much work has been carried out and modern theoretical physics has led us to consider as more 
probable the existence of a multidimensional universe. Recently, the idea was reintroduced in the context of super- 
strings and brane-world theories. This renewed interest can be explained by the fact that multidimensional universes 
may adequately describe known forces and particles and also explain the hierarchy between the gravitational and 
electroweak scales [2,3]. These scenarios have been extended much in recent years and there is now some approaches 
suggesting that many connected parallel branes could exist in an extended bulk (thus creating a so-called "manyfold" ) . 
As summarized in Ref. [4], the existence of such a multi-sheeted spacetime could shed some light on a number of 
L" ■ puzzling cosmological problems including the nature of dark matter structures (which would be identical to normal 
matter but located in a distinct sheet), as well as their invisibility (the gauge fields and most notably electromag- 
netism would be confined within the branes such that the structures belonging to distinct branes would be mutually 
^ ' invisible). 

Following this line of thought, other recent developments have tried to improve the original Kaluza's idea of 
unification by considering discrete extra-dimensions instead of continuous ones. In these models, each discrete point 
in the fifth dimension is endowed with its own metric field such that the resulting universe appears like a multi-sheeted 
spacetime composed by a collection of parallel 4D sheets in interaction. Although such latticized models are recognized 
to suffer from drawbacks (potentially curable), they exhibit promising phenomenological properties, in particular in 
cosmology [5,6]. 

Discrete extra-dimensions have also been studied in the spirit of non commutative geometry. One of the most 
promising approaches is that of Connes and Lott [7], Viet and Wali [8]. Their formulation considers a product 
manifold comprising a continuous four-dimensional part times a discrete one, typically M4 x Zi. The non trivial 
feature of their theory is the coupling of the two spacetime sheets which results in the Higgs field. In several aspects, 
this construction appears as a reminiscence of the five dimensional Kaluza-Klein model with the fifth dimension 
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restricted to only two points. 

In recent papers, present authors have developed a phenomcnological model of a five dimensional two-sheeted space- 
time and studied the quantum behavior of massive particles in such a universe. Different mathematical approaches 
were used involving either a non-commutative product manifold [9] or a two-point discretization of the fifth dimension 
[10]. In the latter approach, the compact extra dimension was treated through a naive discretization scheme involving 
a finite difference analysis. The procedure was inspired from " multigravity" theories where massive gravitons result 
as a consequence of a latticized extra-dimension [5] . Extending such a procedure to the case of the five dimensional 
Dirac equation was demonstrated to provide several original results like fermionic matter oscillations between the 
two four-dimensional sections [9-12]. Therefore, despite all roughness of the approach, the prediction of possible new 
quantum phenomena made the theory interesting to study. 

In the present paper, we propose to extend further the discussion of Refs. [9-12] by considering an asymmetrically 
warped background instead of a flat one. Since each four-dimensional section is now endowed with its own warp 
factor, different physical length scales can be defined on the sheets. We show that under some circumstances, a 
massive particle can oscillate between the two four-dimensional sections and exhibit hyperfast velocities from the 
perspective of a four-dimensional observer thanks to the different length scales in the two sheets. This result is similar 
to that obtained by Chung and Freese in Ref.[13] except that it is here derived by using quantum mechanics exclusively 
without recourse to geodesic calculation. 

The paper is structured as follows. In section II, we derive the two-sheeted Dirac equation for the chosen asymmet- 
rically warped background. The solutions of the Dirac equation are given and compared with the usual ones which 
are recovered at the decoupling limit (infinite distance between the sheets). Then, in the third part, we demonstrate 
the oscillatory behavior in the case of a positive energy particle and the possibility of hyperfast motions. In the fourth 
section, the non relativistic limit of the Dirac equation is derived and the quantum dynamics is studied at low energies. 
It is demonstrated that the oscillatory behavior and the hyperfast velocities still survive in the non relativistic limit, 
thus enabling a possible experimental confirmation of the model. 

II. TWO-SHEETED DIRAC EQUATION IN A WARPED BACKGROUND 

A. Mathematical framework 

In this part, we shall derive the two-sheeted Dirac equation. We start by considering the usual 5D covariant Dirac 
equation: 

(i /D - m) V> = ij A (x, A) p A ij) - rmp = (1) 

with A e {0, 1, 2, 3, 4} and where ^ A {x, A) are the curvature dependent Dirac matrices. Pa is the covariant derivative 
given by : 

Pa = d A + r A (x,X) (2) 

with Ta(x, A) the spin connections and 8a = d/dx A . In Eq. (1) and (2), "x" refers to the four-dimensional coordinates 
and A to the transverse extra dimension (a subscript or an exponent equals to 4 corresponds to A). Note that 
throughout this paper, we will work in units where c = 1, fi = 1. Now, we assume that the five dimensional metric 
takes the following specific form 

ds 2 = g AB dx A dx B = dt 2 - R 2 (X) [dx 2 + dy 2 + dz 2 ] + d\ 2 (3) 

where R (A) is the warp factor. Note that we consider that the warping only involves the space interval of the line 
element and not the temporal part. This is similar to the choice made by Chung and Freese in Ref. [13] and differs 
from the classical brane-world approach where the metric takes a canonical form. Our choice of metric corresponds 
to an asymmetric warping as defined in Ref. [13]. Since the results of this paper are not affected by the exact form of 
R (A), we are keeping this general form throughout this paper. Note that a signature (+, —,—,—,+) is chosen instead 
of the usual (+, — , — , — , — ) one. It is straightforward to show that a timelike extra-dimension is necessary to ensure 
energy conservation in the present model (a spacelike extra-dimension leads to a non-hermitian Hamiltonian as it is 
obvious from Eqs.(lOl) and (102) hereafter). 

The five dimensional Dirac matrices in curved space take the form: 



1 A (x,\)=e A (x,\) 1 a 



(4) 
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where define the vielbein according to: 

g AB = e A (x)e B (x) V ab (5) 

with gAB the 5D metric and r] a b the five dimensional metric tensor of the Minkowski spacetime. The vielbein is given 
by: 

C A< > ../.111 

which leads to 



ef(x) = diag [1, 1) (6) 



'fix, A) = 7°, l l {x, A) = I 7 % 7 4 (x, A) = 7 5 (7) 

where 7 5 = i7°7 1 7 2 7 3 anti-commutes with the usual Dirac matrices 7 M in flat space, such that one verifies {7°, 7 6 } = 
2rj ab . The spin connection must satisfy the expression [l8j |: 

I^z, A) = i 7 s [9a7 B (^, A) + rg A 7 C (x, A)] (8) 

where Tq A are the Christoffel symbols for the metric field defined above. From (8), the following covariant derivative 
terms can be found : 

p = d , A - a*, A = di - (9) 

where a dot implies a derivation along the fifth dimension and the gamma matrices in Eq.(9) are those of flat space. 
The five dimensional Dirac equation for the metric field (3) finally becomes : 

|*7°0o + ^1% + i 7 5 04 + 4| 7 5 - m j V = (10) 

Let us now discretize the fifth dimension. The way to proceed is inspired by those described in Refs. [5-6]. However, 
instead of considering a one dimensional lattice, we focus on a more restricted compact extra-dimension containing 
two points located at coordinates A G {— 5/2, +6/2} (see Ref. [10] where such an approach had been developed). 
At each site, there is thus a four-dimensional submanifold X± where the particle wave function takes the local form 
ip{x,X — +8/2) = 4'±(x). In the proposed geometrical framework the derivative 84 along the discrete dimension 
simply reduces to a finite difference [10] (see Appendix A) : 

0#l ± -» 9 W>± - V> T ) and R ± =g (R± - R T ) (11) 

where g is the inverse of the distance 5 between the four-dimensional sections and R_ = R(X = —5/2), R + = R(X = 
+5/2) are the "projected" warp factors on the two sheets. 

Using those expressions, the five dimensional Dirac equation breaks down into a set of two coupled four-dimensional 
like Dirac equations: 

*7 O 0o^ + + 7^7*0^+ + H 5 g (V>+ - V-) + rfj^9 ( R + - R -) ^+ - = ( 12 ) 

i^do^P- + ^-7*3^- + rfd (V- - iM + *7 5 tt|-.9 (R- - R+) - = (13) 

This equation can be rewritten in a compact matrix form such that 

(i p - m) ip = (14) 

with 

P = T° do + |r" d v + g T 5 + g A (15) 



provided that 



*-i zv-it d-h 1 '* u«-y=( o\°o- 

r ~ 1, -7 5 7 5 J ' 2 I 7 5 - R+) /R- ) • ( ?) 

Eq. (14) is the two-sheeted Dirac equation for the metric field (3) which will be studied in this paper. Note that it 
takes a form which reminds the usual Dirac equation except for the last two terms igT 5 and igA which generate the 
coupling between the sheets. It can be noted that no interaction occurs anymore at the limit where the sheets are 
infinitely separated, i.e. g — ► 0. 

B. Asymmetrical two-sheeted Klein-Gordon equation 

To solve Eq.(14) it is mandatory to introduce the auxiliary field (f> given by : 

V> = P + m)<j) (18) 

= do + % -^r d, + ™ + i g r 5 + % g a 

Then, it can be shown that <j> is solution of: 

(i p - ro) (i P + m) <f> = |-a 2 + ±8$ - m 2 - g fP d t - g 2 r| = (19) 



with: 



H - ^ - 7 4 7 5 (l/i2- - J ' 1 - ^ -n £_ J (2Uj 

where we have set : 

29Rl - 30R+R- + 9R 2 _ 29R 2 - 30R+R- + 9R 2 + 3 (R+ - R-f , . 

s + = 4 ' s = ~ ^ — ' n = 2 2 W (21) 

Following the usual procedure, we can solve this two-sheeted Klein-Gordon equation by using an ansatz of the form 
= 4>Qe~ l ( Ept ~ p ' x \ Introducing this ansatz into Eq.(19) leads to the system: 



E l - 1% - m 2 - ,g 2 S+ g 2 Ii - ig (1/R+ - 1/R_) Wi* 
g 2 U - ig (l/i?_ - 1/R+) 7 VPi E 2 ~j^-m 2 - ,g 2 E_ 



b = (22) 



C. Free-field solutions of the two-sheeted Klein-Gordon and Dirac equations 

Let us determine the solutions of the two-sheeted Klein-Gordon and Dirac equations. To simplify notations in the 
forthcoming results, it is convenient to write: 



R + + R - ^ _ „,2„4 
R 2 + R 2 _ 



X= U»2 and~=aV + /V+7 (23) 



with 



■5 



7 = a 



4n 2 + (£+ - £_ 



(25) 



, 1 1 

« = 9 n, t = z.g — — 



From Eq.(22), the energy eigenvalues E and eigenvectors </>o are easily found: 

1/2 

and 



— — /• with (j)Q = 1 



r}(j>; 



2, iX . 2 ( S + + S - 

"» + P g + - 9 2 



(26) 



(27) 



(28) 



with 77 = | \ap 2 - g 2 (E + - E_) - \/S 



and 



E = E p = Jm z +p*$+g 



X , 2 (E + +E_) VS 



1/2 



r— /• with 0o = 



#A 



and 



with, = (;-) 



(29) 



(30) 



We also have: 



and 



with fj = | ap 2 — g 2 (E + — E_) + \/S . A = ±i and refers throughout the paper to the helicity states of the particle. 



A = k + 2rAp with A = ( X q A ) 



= k — 2rAp with tp\ 





Aa 



where x\ is the spinor such that, if one considers the impulsion p, it verifies the usual equation: 

a l PiX\ = 2Apxa 

One notes that {i)\,<f>\} and <p\} are the eigensolutions of the operator T given by: 

K + TO 1 Pi 



T = k — T7 5 7>i = 



K — TO % Pi 



(31) 



(32) 



(33) 



(34) 



and which constitutes the non-diagonal terms of the matrix in Eq. (22). 
Using Eq.(18), it can be shown that: 



/ E + m -^+igG+ 







-ig \ 



+ igQ + -E + m -ig 





-ig 



with 



-19 




e± = i + 



E + m -?^+igO- 
^+ ig O_ -E + m J 



3 R± — Rzf 
2 R± 



(35) 



(36) 
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For the positive energies, Eqs.(27) and (29) suggest to look for solutions of the form: 



30 = V 



Similarly for negative energies, Eqs.(28) and (30) suggest to look for solutions of the form: 



E>0 = V 





Xx ' 




■ o - 




Xx ' 




■ o ■ 







+ <?A 










+ #A 





V 





Xa 


or (f) Q = fj 





Xx 




. 




. 




. 








(37) 



" o ■ 




- o ■ 




- o ■ 




" o ■ 


XA 




+ r x 






or 0o = fj 


XA 




+ r 






. 




- XA . 




. 




- XA . 



(38) 



Inserting the ansatz (37) and (38) into Eq.(35) leads (after normalization) to the following solutions for the Dirac 
equation: 



For E = ±E V 



ux ip) = -j= 



7j (E p + m) xx 

ri(%f+igQ + )-ig#x) Xx 
i? A (E p + m) xx 

{$ x {^P + igQ_j-i gV j XX 



V X (p) = -7= 



(? (tf ~ lge+ ) + l ^*x)xx 
r)(E p + m) xx 

0i(*£-igQj+igri) XX 

d* x (E p + m) Xx 



(39) 



where u\(p) and v\(p) refer respectively to the positive and negative energy solutions. 



For E = ±E V 



ux (p) = 



fj {E p + m)xx 

fj(^+ig® + ) -ig-&x) \\ 



$x [E p + mj xx 

(#x(%£+igQ-)-igvi \> 



> , vx (p) 



V (^f ~~ + ^a) Xa 

fj (e p + mj xx 

(*a (it-igQ-j+igfi) xx 



(i! p + mj 



Xx 



(40) 



where, as previously, ft\(p) and v\(p) refer respectively to the positive and negative energy solutions. In this model 
the status of the negative energy is the same that in the classical Dirac equation. The normalization constants C and 
C are easily calculated: 



and 



C = V 2 [(E P + mf + |j- + g 2 (1 + e 2 + ) | + \6 X \ 2 [(Ep + mf + |1 + g 2 (l + 6 2 _) j - 2g \9 X \ 
C = V 2 [(e p + m) 2 + +g 2 (l + 6 2 + ) J + \8 X \ 2 [(e p + mj" + |J- + g 2 (l + 6 2 _)| - 2fj\9 x 



(41) 



(42) 



where \tf x \ 2 = k 2 + £ 2 p 2 with £ = g (-±- - -^j . 

It can be easily checked that at the limit where R± — > 1, i.e. when both sheets have identical warp factors, the 
above solutions conform to that of the paper [10]. 

Reciprocally, at the decoupling limit (i.e. when g — > 0), the solutions (39) and (40) become "two times" the 
classical Dirac ones expressed using the coordinates "x". A careful calculation shows that the spinors of positive 
energy become: 



ux = 



y/2E p (Ep + m) 



(E p + m)xx 

2Xp 






and u\ 



2E, 



'p (e p + mj 



(e p + mj xx 



2Xp 
R- 



(43) 
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while the energies become : 



E p = Jm 2 +P 2 -^2 and E p 



'to 2 +P -jjT 



(44) 



These solutions describe two four-dimensional spinors living in distinct sheets provided that the physical momenta 
p + and p_ arc p + = p/R+ and p- = p/R- on the (+) and (— ) sheet respectively. 

This choice implies that the physical length coordinates are x + = R + x and x_ = R-x on the (+) and (— ) sheet 
respectively. Note that under such a coordinate rescaling, both action and phase are conserved since: 



p • x = p+ • (R+x) = p • (R _x) = p+ x+ = p x 



(45) 



More clarifications on the coordinate rescaling will be given later in section IV. Throughout this paper, x and p 
will refer to "global" coordinates and momenta whilst x± and p± will refer to "physical" or "ordinary" coordinates 
and momenta. 

III. FERMIONIC OSCILLATIONS BETWEEN THE SHEETS AND HYPERFAST DISPLACEMENTS 
A. Particle motion between the two four-dimensional sections 



The predictions of the model can be illustrated by studying the following state corresponding to an unpolarizcd 
particle of positive energy: 



V> = (u 1/2 e l P- x e-' E " t - u 1/2 e ip - x e- iL "' 

2\/V ^ 



) + -i= (u_ 1/2 ^ x e- iE ^ - u_y 2 ^ x e 



(46) 



By virtue of the two-sheeted structure of spacetime, tp is a 8-component spinor whose first fourth components are 
located on the (+) sheet and the last fourth on the (— ) sheet. 

The probability P + (resp. P_) to find the particle in the (+) sheet (rcsp. (— ) sheet) is simply given by the 
integration of |V>+| 2 (resp. |?/>-| 2 ), the square norm of the first (resp. last) fourth components of i[>, over the space 
coordinates on the volume V. We get 



P + = A-2B cos (E p - Ej}j t 



and 



with 



P_ = 1 - P \ 



+ g 2 \$\ 2 + 2g V 



+ g 2 \d\ 2 + 2gg - g K Q + 



and 



B = 



1 



2VCC 



'in 



2 

(E p + m) (e p + m) + + g 2 2 + 



+ g 2 \^\ 2 + 2g[7 1 + fj] 



R, 



(47) 



(48) 



(49) 



(50) 



The presence of the cosine term in Eq. (47) shows that the particle oscillates between the two sheets with a frequency 
v proportional to the difference of the two energy eigenvalues, i.e.: 



2-kv = AE = 



E p — E p 



(51) 



If one assumes that the mass term prevails over the other contributions, then the frequency is given by (at first order): 



2tw ^ 



2to 



(52) 
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and it can be easily checked that the oscillation frequency increases with the coupling strength between the two sheets 
and it is enhanced for low mass particles. At the limit of very low impulsion and for identical warp factors (R+ = R-), 
we get v = g 2 / (7rm) which is the result found in Ref. [10] in a flat background. 

Eq.(47) tells us that when the particle oscillates, it can disappear periodically from the perspective of any four- 
dimensional observer. This result is just the discrete counterpart of motions through the bulk predicted in some 
brane- world theories. A common feature of all those models (be they involving continuous or discrete extra-dimension 
[9-12,14,15]) is the apparent violation of the energy conservation from a four-dimensional point of view. It is obvious 
that such a violation is only an artifact of low dimensionality. 



B. Asymmetrical warp factors and velocities 

In order to illustrate the incidence of the different physical length scales on the two sheets (x+ ^ X- in the more 
general case for a same value of x), it is convenient to rewrite the equation (14) by multiplying it by T°. Then, the 
two-sheeted Dirac equation can be recast as 

= Hyj (53) 
where H is the two-sheeted Hamiltonian. In this form, the global velocity operator V can be trivially calculated: 

with p the global momentum operator. The global velocity can then be expressed: 

v = {i/>\V\i/>) = ^-(V + |7Ve i |^+) + ^-^-|7Vei|^-) (55) 

It is instructive to consider the particle velocity measured by physical observers of the two sheets. For an observer of 
the (+) sheet, any particle confined to this sheet moves according to the velocity operator V + (notice the disappearance 
of the warp factor terms): 



oh\ ( 7 <y 



Similarly any particle located in the (— ) sheet will be observed by local observers as moving according to the velocity 
operator: 



Therefore, the global velocity operator of a particle unrestricted in its motion can be written: 
such that in the most general case, the global velocity reads: 



v = v+/R+ +v_/i?_ (59) 

with 

v + = (V|V + |V) = (V + |7Ve i |V+) (60) 

and similarly for v . 

We see that the particle velocity can be expressed as a combination of the physical velocities of the particle, had 
it remained confined in the (+) or the (— ) sheet. As the particle oscillates from one sheet to the other one with a 
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time period T, the effective distance x that the particle is able to travel during one period is given by (assuming a 
rectilinear motion): 



x = 



T 

vdt (61) 



o 



Let us define a starting point x + = (or x_ = 0) at t = in the (+) sheet (or resp. (— ) sheet). Therefore, an 
oscillating particle detected at time t = T in the (+) sheet (or rcsp. (— ) sheet) is observed at a physical distance x+ 
(or resp. x_) from its starting point: 

T T 

x + = R + x = j + ~^T~ V -^J dt and/or x_ = = J + v + ^j dt (62) 

o o 

As a consequence of the oscillatory motion, the particle is observed at time t = T and at a distance x + (or resp. X-) 
as if it had traveled at a fictitious mean velocity v + (or V-) from the perspective of an observer of the (+) sheet (or 
resp. (— ) sheet): 

T 

1 I f 

v+ = -x+ =R + - vdt (63) 
o 



and/or 

1 



T 

1 

x- = R 



f J vdt (64) 



T 

o 



with obviously: v+ ^ v+ and w_ ^ V-. 

Therefore, in the most general case, the fictitious velocity (v±) of an oscillating particle differs from that it would 
have, had it remained confined in only one sheet (i.e. v±). To clarify this issue, let us consider the case of an 
unpolarized particle initially mainly localized in the (+) sheet and described by the wave function 

V> = -^-j= (u 1/2 e lp - x e- lE " t - u 1/2 e lp - x e- 1 ^ + ^= (u^^e 1 ^ ''e~ iE ^ - u_ 1/2 e lp - x e- lE ^ (65) 

Considering that the momentum p is along the Oz axis, and according to the convention used for the spinor deter- 
mination, the particle velocity becomes : 

v = ve 3 = (tp\ V 3 |V>) = J ^V 3 tpd 3 x (66) 

We thus have 

v=^{a-2bcos(E p -E p + <p)t} (67) 



with: 



a = u\ /2 V 3 u 1/2 + iJ_ 1/2 V 3 u_i /2 + u\ /2 V 3 u 1/2 + u J< _ 1/2 V 3 U-i/ 2 (68) 



be tip = u\ /2 V 3 u 1/2 + v)_ 1/2 V 3 u_ 1/2 and be ' Llf = u\ /2 V 3 u 1/2 + u j *_ 1/2 V 3 u_ 1/2 (69) 
As a consequence, the mean fictitious velocity as seen by an observer of the (+) sheet becomes: 

T 

v+ = J vdt = - A R+a (70) 

o 
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After tedious calculations, the fictitious particle velocity from the perspective of an observer located in the (+) sheet 
reads: 

v+ = (E p + m) 

For illustrative purposes, it is convenient to simplify this expression by considering its non relativistic limit. Then 
assuming that g « m, p and p « m, one gets the more compact form 

"^{^ + °[») 2 }" +o w 2 =K i+ S)^ < 72 > 

with x given by Eq.(23). 

On the other hand, at the decoupling limit g — > and assuming completely localized particles in the (+) sheet the 
velocity can be calculated to be (from Eq.(44)) 



R\ 



K/2 

Rl 



R+^L ( E p 
C V 



+ m 



|^l/2 

R?_ 



m 2 



(71) 



v+ = — — = — (73) 
E p R+ E p 

At low impulsion, this expression conforms with the usual one as expected, i.e. v + = p + /m. As a consequence, the 
expression (72) can be conveniently rewritten as 

B +=K i+ t)- + (74) 

We see confirmed that an oscillating particle can travel between two locations with a fictitious (but effective) speed 
v + which differs (in the general case where R+ ^ R-) from that it would have, had it remained in only one sheet (i.e. 
v + in the (+) sheet). If J?_ << R + , the apparent particle velocity v + can become huge even if the particle velocity 
v + in the (+) sheet remains moderated. 

A numerical example can be given for illustrative purposes. Let us assume a warp factor ratio = 250 and 

a particle initially located in the (+) sheet with a non relativistic velocity v + = 20 km ■ s _1 . Note that since we have 
v+/v- = R_/R + in the non relativistic limit (from Eq.(73)), we see immediately that if a particle initially in the 
(+) sheet reaches the (— ) sheet and stays there, its physical velocity as measured by an observed of the (— ) sheet 
will not be 20 km ■ s _1 but 5000 km ■ s _1 . Although this is an important velocity increase, this value is still a non 
relativistic one. At this stage, it is important to stress that the equations of the physical velocities t> + (see Eq.(73)) 
and V- imply that v± < 1 whatever R + , p, m. Therefore, even if the particle moves faster in the (— ) sheet, it 
can never exceed the light velocity. The laws of special relativity are safe in the present approach. In addition, from 
Eq.(45) we note that x+/x- = R+/R-. In the present example, this means that the distances are 250 times shorter 
in the (— ) sheet than in the (+) sheet. Let us consider a particle initially localized in the (+) sheet and transferred 
into the (— ) sheet. The particle covers a distance during a time t and goes back in the (+) sheet. For an observer 
of the (+) sheet, the detected particle has apparently moved a distance x + during a time t, such that its apparent 
(but fictitious) velocity v+ = x+/t or v+ = v+ (R+/R-) ■ The latter expression is the consequence of the shortened 
distance and increased proper velocity in (— ) sheet. The reason why this expression differs from Eq.(74) is the result 
of the oscillating behavior of the particle. Indeed, on the average, the oscillating particle spends as much time in sheet 
(— ) than in sheet (+). The fictitious velocity v + in Eq.(74) is then an average velocity between the real velocity in 
(+) sheet and the fictitious velocity, had it remained in the (— ) sheet only. 

For a proper velocity v + of 20 km ■ s _1 in the (+) sheet, Eq.(74) shows that the fictitious velocity v + (calculated by 
an observer located in the (+) sheet) is about 2 times the light velocity. Although the proper velocity of the particle 
remains non relativistic in both sheets, its apparent velocity from the perspective of an observer of the (+) sheet 
exceeds now the light speed. 

In Ref. [13] , it was suggested that the homogeneity of the universe by the time of nucleosynthesis might be explained 
by particles motions in an asymmetric two-branes system (one of which being a "hidden" brane). According to the 
authors, an impulse originating on one brane of the system can take a shortcut through the other brane and affect 
our brane at a point outside the conventional causal horizon. The present paper tries to go further on those aspects 
by using a very simple although realistic quantum mechanical model. Although our approach is radically different 
from that of Ref. [13], our results share obvious similarities (e.g. the particles confined in the other sheet are invisible 
to us, the oscillating particles can reach distant point outside their "naive" horizon...). Hence, it is believed that 
there may generically exist non-inflationary solution to the horizon problem in theories with extra-dimensions (be 
they continuous or discontinuous). 
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IV. TWO-SHEETED ELECTROMAGNETIC FIELD AND PAULI EQUATION 
A. Introduction of the Electromagnetic field 

To account for the discrete structure of the bulk, the usual U(l) gauge field must be substituted by an extended 
U(l) ® U(l) gauge field such that 

q _ l4x4exp(-zgA + ) 

14x4 exp(— iqA-) 

We look for an appropriate gauge such that p A — > p+ /A with the following rule of transformation 

/A' = G + G[p,Gft] (76) 

A convenient choice is (see Refs. [9-12,16]) 

/A= 7 5 x 

' 7 5 x f njc^A- 

where 7 M and j 5 are the usual Dirac matrices. Choosing x = whatever the gauge choice can be done only if 
A+ = A_ = A. 

In addition, the above definitions impose the following gauge transformations 

A± = A± + d^A (78) 

On the two sheets live the distinct A + and A_ fields. Each spacetime sheet possesses its own current and charge 
density distribution as sources of the local electromagnetic fields. The off-diagonal term has been set equals to zero. 
This free choice allows a further simplification of the model. If this term is different from zero, it leads to a coupling 
between the two photons fields such that each charged particle becomes sensitive to the electromagnetic fields of both 
sheets irrespective of its localization in the bulk. With the present choice, the electromagnetic field of a sheet couples 
only with the particles belonging to the same sheet. 



(75) 



(77) 



B. Derivation of the Pauli equation 



Let us derive the non relativistic limit of the two-sheeted Dirac equation for the metric field (3) . We first introduce 
the gauge contributions /A of both sheets into the two-sheeted Dirac equation following the standard procedure 



such that 



One can also write 



where 



(ip A - to)* = 
Pa = r° (do + iqAo) + (d v + iqA^j + gT 5 + gA 
ido^ = -ir°-^r" (dr, + iqA^ * - igT°r 5 y - igT°A^ + mT ^ + qA f 





A~ 



(79) 
(80) 

(81) 

(82) 



When m is large compared with the kinetic energy, the most rapid time dependence is in the factor exp(±zmi). For 
a free positive energy particle, and for small kinetic and electromagnetic energies, we may therefore seek a solution of 
the form * = %pe~' lmt with 



X- 
9- 



(83) 
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where x+, 9 + , X-, 0- are two-component spinors. One can then write 

id oX + = qA+ X + - i(l/R+)<r v ( d v + °+ ~ ^ ( + ~ - W (3/2) {(R+ - R-) /R+} 0+ (84) 

idox- = qA oX - - i(l/R-)<r v {d v + iqA~) 0- + ig (0+ - 6.) - ig (3/2) {(R_ - R+) /R-}9- (85) 

id o 0+ = -i(l/R+)a v (8 V + iqA+) X + + qA+9+ + ig ( X+ - X-) - 2m9 + + ig (3/2) {(R + - R_) /R + } X+ (86) 

id o 0- = -i(l/R-)a v (d n + iqA~) X - + qA 9- - ig ( X + - X-) - 2m#_ + ig (3/2) {(R_ - R+) /R-} X - (87) 

In addition, when p is much smaller than to, + and 9_ become tiny in comparison with x+ an d X-- For small 
electromagnetic and kinetic energies, one then get the following expressions 

9+ « -i^/R^l-ar, (8 n + zqA+) x+ + (X+ - X-) + ^ (3/2) {(R+ - R-) /R+} X+ (88) 

« (9, + zgA-) x- - (X+ - X-) + (3/2) - R+) /R-} X- (89) 

In order to give these equations a more conventional form, vectors must be used instead of covariant or contravariant 
terms. The procedure, which is similar to that used in cosmology [17], can be summarized as follows. For a metric 
field given by gij = R 2 5ij = 1,2,3), we can define an orthonormal basis: 

Bi = | = Rg (90) 

such that g^ = gj • gj = g^. What are usually called the components of a vector a in elementary treatments arc 
neither the covariant components at nor the contravariant components a" 1 , but the "ordinary" components: 

at = a • = Ra l = R^ 1 ^ (91) 

Moreover, to be consistent, vectors and operators relative to a specific sheet have to be expressed using the metric of 
the corresponding sheet, i.e. following Ref. [17] 

v ± = E e ^^7 = ^Ep* = *£p^ where p* = and p±.< = -^^r = ^ (92) 

and similarly, magnetic vector potentials and magnetic fields are given by 

A±,i = R±Ai_ = R± 1 A± ii and similarly for B± d (93) 

Note that we have introduced the "global" vector p exactly as we did for the Dirac equation (Eq.(45)). 

By setting $± and A± for the usual electric and magnetic potential, it can be easily shown that the Pauli equation 
reads 



Hip = ido<p with ip = 

where we have used the fact that 



X+ 
X- 



(94) 



a 7] a v {V± - iqA±) r} {V± - iqA±) I/ = (V± - iqA±) +qa-B ± (95) 
The resulting Hamiltonian can be written as the following sum 

H=(H k + H m + H p )+(H c + H cm + H cp ) (96) 
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where 



Hk = ~- — 



1 

2m 



(V+ - iqA + Y 









(V_ -iqA- f 



H„ 



q 

2m 



a ■ B+ 
<t-B 



(97) 
(98) 



Hp = 



q$+ 
q<P- 



(99) 



H c = 9 - 

m 



1 + (3/2) (2*0=) + (9/8) (2*0=)' 
_1 + (3/4)1^1! 



1 + 



(3/2) (> 



R 



(3/4) 

^) + (9/8) (I 



-R+ 



(100) 



2m 



cr-{A+-A_} 
-cr-{A + -A_} 



(101) 



TT _ 

cp " 2m 



cr-{V+-V_} 
-ct-{V+-V_} 



(102) 



This two-sheeted Pauli equation is very similar to that derived in Refs [9] and [10] except for the gradient operator 
which is distinct in the two sheets as a consequence of the two metric fields. It is worth being noticed that at the 
limit of decoupling, we get two times the usual Pauli equation as expected. Therefore, at the limit of small g, the 
difference between the two-sheeted Pauli equation and the usual one are not expected to be significant. 

The first three terms of the Hamiltonian correspond to the classical contributions of the "one-sheeted" Pauli's 
equation. H^ is the kinetic term whereas H m and H p relate to the magnetic and coulomb terms respectively. The 
last three terms correspond to new predictions of the model. The term H c behaves as a constant coupling between 
the two sheets. This term is responsible for the spontaneous particle oscillations studied previously in the relativistic 
limit. H cp and H cm introduce the geometrical coupling between the sheets through kinetic and magnetic terms. It is 
worth stressing that the coupling which involves the kinetic mixing disappears if the two sheets have the same warp 
factors whereas H cm does not. 



C. Hyperfast velocities in the non relativistic limit 

As previously, three different velocity operators can be defined : V + (V_ ) which corresponds to the velocity of a 
confined particle in the (+) (resp.(— )) sheet and the global velocity operator V. These operators read 



dH 



(p + -qA + ) £0 



2m 



2m 




and 



and 



V = 



dH 
0p~ 



2m m 



qA_) 



R+ dp+ 



(103) 



(104) 




(105) 



We point out the existence of a non diagonal components of the velocity operator. However, since this term is 
proportional to g we can ignore it for the moment and in the forthcoming calculations to concentrate only on the 
effect of the diagonal terms in V. 



14 



D. Illustrative case 



In the following, we neglect the terms g 2 of the hamiltonian in front of the g terms. In this way we underline 
the contribution of H cm and H cp . To illustrate the survivance of hyperfast displacement in the non relativistic 
limit, it is convenient to consider the case of a particle initially located in the first sheet and embedded in a region 
of constant curlless magnetic vector potential. For simplicity reasons, let us consider a neutron-like particle, i.e. 
a chargeless particle with a magnetic moment. We do not consider the complications arising from the anomalous 
magnetic moment of this particle. Instead, we assume that it is identical to that of an electron. In absence of any 
magnetic field nor scalar potential and by neglecting the g 2 terms in front of the g terms, the Hamiltonian reads 



H 

with a = g/(2m) and 



K + —iota ■ P 
iaa ■ P K_ 



(106) 



= < 107 > 

and 

p = eA+ {^-_L}p (108) 

In the following, one considers that A and p are colinear. The eigenvectors of the above hamiltonian are 

(109) 



1 



/N± 

with the corresponding eigenvalues 



(E±-K-) X \ 
iaa ■ Pxa 



E± = \ (k+ +K_± sj{K + -K_f + 4a 2 P 2 ) 



(110) 



Xa is a spinor where A = ±1/2 and stands for both spin states and N± = (E± — K_) (2E± — K + — K_). 

It is not difficult to convince oneself that the particle ability of reaching the second spacetime sheet still survives 
in the non-relativistic limit. The wave function assuming that the particle is initially (t — 0) located in the (+) sheet 
with some polarization state {n\^ 2 — n2 _i/ 2 )/( n i/2 + n -i/2) ^ s 

1>=^= N/2 (au +jl/2 e- iJ5 +* + &u_, 1/2 e- iJ5 -*) + n_ 1/2 (an +j _ 1/2 e- iB +* + &u_,_ 1/2 e- ii5 -') } e** (111) 

with n\^ 2 + n 1 ^^ = 1 and 



N+ / N 

±_and6=W- — -— (112) 



y N+ + N- y N+ + N- 

The probability P to find the particle in the second sheet can be trivially calculated 

sin 2 {(l/2)(i; + -£L)t} (113) 



1 -,2 



1 + K 

with 



K+ - K_ 
2aP 



(114) 



Without going further, it can already be noticed that if A = 



15 



and 



P 2 + 

E+ — E— — — — 
+ 2m 



R+ 




(116) 



with p + = p/R+- In absence of confining effect, the particle oscillates between the two four-dimensional sections 
provided that p+ ^ 0. However, if p + becomes larger than g, P quickly drops to zero and the spontaneous oscillations 
are strongly suppressed. 

Returning back to the general case where A ^ and p + ^ 0, it is interesting to calculate the particle velocity. 
Assuming p and A are oriented along the Oz axis, the global particle velocity v is given by 



v = v z = J v>V 2 



ipdV (117) 



1 v I v ^(K+-K_f + 4a 2 P 2 
= ^-ilcos 2 {(l/2)(£ + -£_)i} + - r - -± = sin 2 {(1/2) (£+-£_)£} 

Since the particle oscillates between the two sheets, the velocity exhibits also an oscillatory behavior. The fictitious 
velocity v + from the point of view of an observer in the sheet (+), after a period T = ^ E ^Z E ) is: 



and therefore: 



v = tl± I v dt (116 
T Jo 




From Eqs.(107) and (108) and for a large enough A, the previous expression simplifies further: 



with v+ = p+/m the usual particle velocity had it remained confined in the (+) sheet. We recover the relation already 
derived from the relativistic approach at the limit of low velocities. 

Again, we see that if i?+ >> J?_, the particle moves with an apparent velocity that can exceed the velocity v + 
this particle would have, had it remained in the (+) sheet. The novelty in comparison with the result (74) is that 
the particle oscillates through the application of a magnetic potential and therefore this result suggest that the model 
could perhaps be experimentally investigated (this possibility and the properties of H cm have been explored more in 
Rcfs. [11,12]). 

We stress that the results of this paper have been obtained for a free particle, e.g. by assuming that the magnetic 
field, the scalar potential and any other environmental contribution can be neglected (excepted for a constant curllcss 
magnetic vector potential as discussed previously). This condition is obviously a very restrictive one and in most 
cases, the suppression of these terms could be hardly achieved. In Ref. [12], it was demonstrated that any diagonal 
term in the Hamiltonian strongly suppress the particle oscillations. In fact, it is not difficult to demonstrate that 
the more energetic the particle is, the more the oscillations are suppressed. Even the gravitational potential, whose 
contribution modifies the particle energy, will affect the oscillations and restricts the particle motion between the 
sheets. As demonstrated in Ref. [12], our model suggests that any massive particle could be spontaneously confined 
within the sheets (through environmental interactions) without requiring any complementary scalar field or repulsive 
gravity [12]. Since no particle disappearance has been noted to date, it is very likely that the degree of confinement is 
very strong, and/or that the coupling constant g is very small. Therefore, it is expected that the phenomena described 
in this paper will be hardly observed, especially if one takes into account the fact that no current experimental setup is 
suitably designed for searching for these phenomena. Any experiment aiming at demonstrated the behavior predicted 
in the present paper will require very particular conditions although they might not be completely out of reach of our 
present technology [11,12]. 
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V. CONCLUSIONS 



In this paper, we have studied the quantum dynamics of spin half particles in an asymmetrically two-sheeted 
spacetimc. It was shown that any free particle oscillates between the two sheets as a consequence of the geometrical 
coupling along the discrete extra-dimension. By oscillating, any massive particle is able to travel between distant 
points which are normally outside its four-dimensional horizon. The reason arises from the differential warping which 
leads to reduced length scales and increased velocities in one of the sheet. 



APPENDIX A 



The discrete derivative used in present and previous papers [10-12] is defined as follows. 

A compact oriented discretized dimension Z n — {si\i € {0, n — 1}} is considered, with Sj are each of the n sites 
of Z n (see figure 1). Z n is invariant under the cyclic group Z/nZ. Z n is oriented positively with increasing values 
of site index i. For each site Sj of Z n it is possible to define a coordinate A(sj) = Aj = Si + q. Note that both 5 
and q have the dimension of a length, q is an arbitrary constant (the presence of this term will be clarified shortly 
afterwards) whereas 6 is the distance between nearest-neighbor sites. 




FIG. 1: General polygonal representation of the compact discretized Z n dimension (left) and of its Z-z limit (right). 



One then defines a positive restricted algebraic distance d a (si, Sj) on Z n . That means that one goes from site Sj to 
site Sj by imposing a positive direction, i.e. 

d ^ s o) = { \ . - \ t + 5n if i < \ (A1) 

This imposes 

d (si_i,Si) = S for i e {l,...,n- 1} (A2) 

and 

d o (s n _i,s ) = S (A3) 

with respect to the cyclic properties. 

If one considers an arbitrary function (f)(X(si)) — <p{\i) — (j>i, the derivative d\<fr(\i) at each site Si is naturally given 
by: 

<m(ao ee for . e {1) n _ 1} (A4) 

a a {Si-i, St) 

with of course 

d^Xo) ee ^°;-^ (A "7 l) (A5) 

These expressions can be written in the usual form 

dxfa ee (1/5) {fa - fa-i) for i e {1, ...,n - 1} (A6) 



17 

and 

6> A 0o = {1/5) (0 O - 0„-l) (A7) 

By construction, the derivative is invariant through cyclic permutation and corresponds to the so-called "Z n - 
derivative" . 

The "^-derivative" is simply obtained by setting n = 2 in previous expressions. The result reads: 

d x 0i = (1/6) (0! - 4> ) and ch> = (1/S) (4> - 4>i) (A8) 
Considering the following substitution (0 O , 4>i) — > 0+) the derivative can be expressed into the form 

A 0± = ±(V«) (*+ - (A9) 
In addition, by setting q = —5/2, wet get A± = ±5/2 for the coordinate of each site. 
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